function sys=plant(t,x,U,D)
sys=zeros(2,1);
J=10;

sys(1)=x(2);  
sys(2)=1/J*(U+D);
end